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Abstract. We introduce a new complex formalism to de- 
scribe arclet fields in clusters of galaxies and derive the 
appropriate inversion techniques to find the related mass 
distribution of the lensing cluster. Applying the complex 
formalism to the statistics of lensed background sources 
shows that the mean of the axial ratios of locally lensed 
isotropically distributed elliptical sources is not influenced 
by the intrinsic ellipticity distribution, if the cluster is sub- 
critical in that region. 
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1. Introduction 

A background population of galaxies seen through a clus- 
ter shows a deformation of the morphology of each indi- 
vidual galaxy as imposed by gravitational lensing. If these 
galaxies are located near the caustics of the mapping, de- 
scribing the lensing properties of the cluster, we see very 
elongated giant arcs. However all images of all galaxies, 
even far from the caustics, are at least slightly distorted, 
which are then called arclets. This distortion could be 
used to trace the potential of the cluster. If we assume 
the unlensed galaxies to have small, circular images, grav- 
itational lensing stretches them to small ellipses. The axial 
ratio of these ellipses and their orientation could be used to 
find the potential, and related through Poisson's equation, 
the projected density distribution. In the case of infinitely 
small images the axial ratio is given by the ratio of the 
eigenvalues of the Jacobian of the lens mapping and the 
orientation is found by rotating the Jacobian to the diag- 
onal form. The task is to reconstruct the potential from a 
measured field of axial ratios and orientations. Some au- 
thors try to fit a mass-model to the cluster to reconstruct 
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the observed arclet field, others derive a differential equa- 
tion for the surface density distribution, which describes 
the effect approximately. 

The assumption of circular sources is by no means true. 
The intrinsic distribution of ellipticities and orientations 
are not known and could only by compared statistically 
to fields of background galaxies which are hopefully not 
lensed. It must therefore be clarified in which sense a mean 
source could be accepted as circular in order to apply the 
techniques mentioned above. 

Technically, the authors use cartesian coordinates or 
a vector formalism to solve the problem. We show that 
using a complex approach leads at least to a very elegant 
and compact, exact representation of the problem. 

For an introduction to the theory and application 
of arclets see Miralda-Escude (1991), Kaiser k Squires 

(1993) , Breimer (1994), Broadhurst, Taylor k Peacock 

(1994) , Carlberg, Yee k Ellington (1994), Fahlman et 
al. (1994), Kneib et al. (1994), Miralda-Escude k Babul 
(1994), Schramm k Kayser (1994), Schneider k Seitz 
(1994), Smail, Ellis k Fitchett (1994a,1994b) and recently, 
containing most references and the basic concepts, the re- 
view of Fort k Mellier (1994). 

2. Mathematical preliminaries 
2.1. Gravitational lensing 

The complex gravitational lens mapping is given by 



w(x, y) = (x + it/) + 



d_ 
dx 



(i) 



and, relating the two dimensional gravitational potential 
to the surface density a, 
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that we can express these equations in terms of Wirtinger 
derivatives (see also section "The radial case" below and 
compare Henrici, 1986) ) 
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The lens mapping is generally not analytic, which means 
that the Wirtinger derivative with respect to ~z is not zero. 
So, if we write an arbitrary complex valued function as 
f(z) it should be implied that it could depend on ~z as 
well. We find 
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where we introduced the complex deflection angle a for 
convenience. With Eq. (5) and Eq. (8) we find 
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And for the derivative of Eq. 5 we find using Eq. 7 

dw a 
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2.2. The Beltrami equation 

For any complex valued function f(z) we find at any po- 
sition z an infinitesimal ellipse which is mapped onto a 
circle by /. The ellipse is characterized by the complex 
dilatation n, which is also called Beltrami parameter. The 
axial ratio of the ellipse is given by 



and the argument of the major axis 
2ip = it + arg(//) . 



(12) 
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In Fig. 1 we show in the complex //-plane the correspond- 
ing ellipses. 

The Beltrami parameter itself is determined in terms 
of Wirtinger derivatives of / by the Beltrami : equation (for 
details see Appendix A) 



(14) 



Fig. 1. The ellipse field over the complex (p = m + ira)-plane. 
The unit-circle is shown for convenience. 



Clearly n(z) is a complex, generally nonanalytic function 
over the complex plane. 1 Since the Jacobian of / is 
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dz 



df(z) 



(15) 



J = implies = 1. For gravitational lensing it is inter- 
esting to note that in the matter free case (a = 0) = 1 
and so [i is completely determined by An integration, 
interpreted as anti derivative, lets us recover w up to a 
function depending on z, which is just z if compared with 
the structure of the lens mapping. 

A mapping is said to be k-quasiconformal (q.c.) if 
I// 1 < k < 1 which means that there is a fixed bound on 
the stretching for the mapping in any given direction com- 
pared to any other direction. Lens mappings are therefore 
q.c. if there are no critical lines. It has been shown that 
there is always a unique solution for the Beltrami equation 
if the mapping is normalized at infinity and q.c. . Unfor- 
tunately, this applies only for marginal clusters and the 
proof of existence and uniqueness is hardly usable for the 
numerical treatment (see e.g. Carleson & Gamelin 1993, 
Ahlfors 1966, Henrici 1986, Imayoshi & Taniguchi 1992). 
In our section "Numerics" we show how to find solutions 
for arbitrary lens mappings. We leave it to the mathemati- 
cians (who are invited) to see if our numerical solution 
applies to a larger class of e. g. Lagrange mappings if ex- 
pressed in terms of integral operators. Even existence and 
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information we could get from arclet fields at all. 

2.3. Stokes's and Pompeiu's formulas for lens mappings 

Pompeiu's formula is the generalization of Cauchy's for- 
mula for nonanalytical functions. Sometimes it is also 
called the Cauchy-Green formula, because it can be de- 
rived using Stoke 's formula, which reads using Wirtinger 
derivatives 
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for an area R with boundary dR. If applied to a function 
one finds Pompeiu's formula: 
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The value of a function at a location inside a region R 
bounded by dR is given by a contour integral over the 
boundary plus a surface integral over the region. If we 
assume that 
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which states that the lens acts as a point mass outside of 
the region R, Pompeiu's formula gives by residues for the 
lens mapping 



w(z) = (1 - l/R 2 )z 
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where the l/R 2 term does not contribute, if in the contour 
integral = R 2 — ► oo. 

We add that if Stoke's theorem Eq. (17) is applied to 
the lens mapping in the form Eq. (6) we find using Eq. (11) 
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which could be used for measuring the mass M inside a 
region R if the deflection angle is known on the boundary. 

2.4. The radial case 

We briefly recall how the definition of the Wirtinger 
derivatives is founded. Remember 
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can be written as 

df = — dx + — dy 
dx dy 
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If we insert the total differentials for da;, dy as imposed by 
Eqs. (22-25), we find 
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which leads to the definition given in Eq. (4). We apply 
this to the polar description of complex numbers 
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The total derivative is 
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A (lens) mapping with radial symmetry would read 

w(r, ip) = R(r)e llp . (38) 

For this mapping the polar Wirtinger derivatives read sim- 
ply 
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It is convenient to set 

n = W 2 ^ , 

so we find for the Beltrami equation 
dR R , , fdR R^ 
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find an ordinary differential equation 

dR 1 R _ 
dr e{r) r 

which is easily solved 
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Measuring the variation of the axial ratio leads to the ra- 
dial lens equation. With Eq. (10) and the radial Beltrami 
equation (42) we find for the density 



a = 2 - 1 



1\ R 



(45) 



These results have been found using the normal vector for- 
malism by Schramm & Kayser (1994) (see also Schneider 
k Seitz, 1994). 



2.5. Some instructive examples 
2.5.1. Point mass 

Let us take a pure point mass acting as a lens. The lens 
mapping reads 



1 / r 

w(z) = z — — = r ) e 1 ^ 

z \ r . 



The ellipse field is given by the Beltrami parameter 
1 1 



^ " z 2 ~ r 2 



and the corresponding axial ratios by 
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The ellipse field is given in Fig. 2. The corresponding vari- 
ation of and e, respectively, is shown in Fig. 3. 

2.5.2. Point mass and shear 

We loose the radial symmetry if a shear is introduced 



1 

w(z) = Z + JZ — — 



(49) 



The corresponding ellipse field is given by the Beltrami 
parameter 



and the corresponding axial ratios 

= \z\ 2 -\l + 7 z 2 \ 
' \z\ 2 + \l + 1 z 2 \ • 

The ellipse field is given in Fig. 4 
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Fig. 2. The ellipse field due to a point mass. Note that the 
large principal axes are fixed. For convenience we also show 
the critical curve. 



2.5.3. Isothermal sphere 

An isothermal sphere is given by a density distribution 
which varies as 



a = <7 - 



1 



Our lens equation then reads 



w(z) = Z — <7q- 



A 



and for the complex dilatation we find 
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which could, of course, be derived with the radial approach 
given above. We show the ellipse field in Fig. 5 and the 
related variation of e and in Fig. 6. 

3. Numerics 

Our intention is to introduce the Beltrami approach to 
cluster lensing and not to discuss detailed applications. 
However, we briefly touch some possibilities how to solve 
the Beltrami equation for lensing problems. Applications 
will be published in a subsequent paper. 

The Wirtinger derivatives are well defined as deriva- 
tives in the sense of distributions. However, it is hardly 
possible to give them a constructive numerical sense with- 
out casting them into their components. 

As typical in potential theory we have isolated or 
boundary value problems. Even if we are able to define 



Fig. 3. Norm of fi and e for the point mass lens. The critical 
line is denoted by a circle on the curve for fi. The related axial 
ratio is zero at that point. 



a homogeneous //-field over a region (ccd frame) for the 
latter case we have to guess the values of the mapping on 
the boundary to find a unique solution of the Beltrami 
equation. For this purpose we decompose the Beltrami 
equation into terms containing normal partial derivatives. 
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Introducing the complex function 
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the equation reads quite simply 
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This equation can be descretised and solved as a bound- 
ary value problem involving all the problems of numerical 
solutions of first order differential equations. We do not 
follow this track here but we assume that we are so lucky 
to have all the contributing material in our region with 
defined //-field. So we can treat our problem as isolated. 
This includes the case where the contributions of the outer 
parts cancel, as it is the case for radial or homoeoidal (see 
e.g. Schramm 1994 for details ) mass distributions: If we 
apply the Beltrami equation to the lens mapping Eq. (9) 
we find: 
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Fig. 4. The ellipse field due to a point mass and an additional 
shear. The critical curve is also shown. 



If we now assume a grid of size (n x m) representing an 
area (x2 — x\) x (t/2 — Hi) in the lens plane with pixel size 
dx = (x'2 — x\)jn and dy = (t/2 — y i ) / m . We find the 
center coordinate = (xi +ij/j) of each pixel by 
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2/1 + jdy 



dx 

T 
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We introduce the numbering scheme 
k = (j — l)n + i 
with the inversion 
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For convenience we define a distance (matrix) function D 
on the grid 
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which applies to the Beltrami equation if the integral is 
approximated by a sum over all pixels with constant den- 
sity. We find for each pixel k: 



dxdy -^—^ ui 
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and 



which , after introduction of vectors <t = (ci . . . a„ 
) can be written as matrix equation: 

Do- = n (66) 



Fig. 5. Ellipse field of the lens as given by an isothermal sphere 
with (To = 1. Note the tangential and radial stretching of the 
ellipses near the critical line and the center. 



Fig. 6. The norm of the complex dilatation fi and the corre- 
sponding axial ratio e. Note the pole of \fi\ where the ellipses 
turn from radial to tangential stretching. 



with obvious solution 
<r=D - V • 



(67) 



Note that adding a sheet of constant density k does 
not alter the //-field if every mass element is scaled by a 
factor 1 + k (magnification transformation). However, if 
we are able to observe arclets out to regions where we 
are sure that the local density of the cluster is zero, we 
can uniquely determine the total mass up to the mass 
contribution of the global cosmological background. 

4. Seeing 

To see the influence of the atmospheric seeing we inspect 
again our definition of the Beltrami parameter Eq. (12). 
If we assume the axis of the arclet to be a, & we would find 
without seeing 



a + b 
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With seeing s the ob 
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respectively, if we assume a gaus- 



sian surface brightness distribution of the source. We find 
therefore 



\fb* 



(69) 



We see that we can recover the //-parameter if the seeing 
is known. However, we see also that we run into problems 
if the seeing is very bad, since the expression for /j, gets 
indefinite. 



5. Statistics 

Observed arclets are, as mentioned above, in general not 
images of circular sources. Additionally, the sources are 
distributed over a redshift range. Although the redshift 
distribution can be accounted for in the Beltrami equation 
we here assume for simplicity that all sources are at the 
same redshift. 

However, we can also describe the situation including 
intrinsic ellipticities using the Beltrami equation if a com- 
posed mapping h is introduced, h is a composition of two 
mappings w and u>. The mapping w maps the observed 
elliptical arclets to the unknown real elliptical source and 
ijj maps this source to a circle. (See Fig. 7). 

The corresponding Beltrami parameters are // w and \i u 
respectively. 



h(z) = w(w(z)) 

which all depend on complex conjugated variables as 
we find therefore for the total differential 
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For the Wirtinger derivatives we find 
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Fig. 7. The composed mapping h maps observed elliptical ar- 
clets to circles. It is composed from the lens mapping w that 
maps the observed arclet to the real elliptical source and an 
artificial to mapping that maps the real elliptical source to a 
circle. 



and so for the Beltrami parameter 
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For simplicity we use indices in parentheses to denote 
derivatives. The last equation then reads 



(77) 



If we suppose now that we are in the linear regime of the 
mapping, we can always write 



Ul = W + /Uu.W, < \fl w I < 1 

w = (1 — <j)z + *fz 
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Note that we do not have to restrict the lens equation to 
the case of real j. However, the only difference would be 
a constant phase in the appropriate Beltrami parameter, 
which could also be attached afterwards. 
Applying Eq.(77) we find 
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and the inverse 
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These equations describe the transformation of the Bel- 
trami parameter under local lensing, given by ;U m . We are 
now interested to see how a local distribution $ of 
transforms under local lensing. This means how the sta- 
tistical properties of a field of elliptical sources alter if 
mapped by a lens locally described by constant density 
and shear: 
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where we have to insert for J w the Jacobian of the mapping 
for \i u 
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Assuming now isotropy for the original sources, we find 



1 
2w 



(1 + /Ift/lw) 



^(MMI) 



(91) 



The normalization for polar coordinates reads as fol- 
lows: 
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which is for uncorrelated and <p 
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For an isotropical distribution 
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and for the radial part 
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t*w = l/Uu.le 1 ^'* 



which is in case of isotropy just 
(t* u ) = 
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because of the angle integration. The image /«h((/«ui}) is 
therefore 



1 + A«w (/"a.) 

as expected. The mean of //^ is 
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This yields for isotropy 
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The integral in square brackets could be evaluated if cast 
into real and imaginary part. The important fact to note is 
that its result depends on sign ((/«wA«ui — !)/(! + A*wA«ui))- 
We use the theory of residues and cast the integral into a 
contour integral around the origin: 
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We therefore find for [jih) using the normalization Eq. (95) 
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Since l/^l < 1 the last condition is fulfilled for |// w | < 1, 
which is valid for a region where the lens mapping is ori- 
entation preserving (quasi conformal). The bad news is 



mation, the mean eccentricity of the observed arclets is in- 
fluenced by the distribution of the intrinsic eccentricities. 
The good news is that we can conclude the intrinsic dis- 
tribution from regions where the mapping is surely quasi 
conformal. 
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A. Appendix A: Why complex? 

The complex notation is typically only used for mappings 
which are analytic. However, even for nonanalytical func- 
tions there could be a remaining structure, which makes 
complex functions easier to use than the 2D-vector ana- 
logue. This has e.g. been shown for the calculation of 
caustics (Witt 1990) and for the computation of elliptical 
lenses (Kassiola & Kovner 1993, Schramm 1994). We show 
now why this is also the case for cluster lensing (compare 
Ahlfors 1966, p4). 

Suppose an arbitrary local (linear) map given by an 
expansion of w = w(z), with w = u + iv and z = x + iy 



v (x) x - 



■ u (y)y 
v (y) y 



(Al) 
(A2) 



where, as above, subscripts in parentheses denote differ- 
entiation. A local ellipse is then mapped onto a circle as 
given by the squared deviation function SDF (compare 
Schramm & Kayser, 1987) 

u 2 + v 2 = (u( x) x + U( y) y) 2 + (v( x) x + V( y) y) 2 = const. (A3) 
which reads in classical notation 
u 2 + v 2 = Ex 2 + 2Fxy + Gy 2 = const, 
with 
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1 (X) T V (x) ! 

l( x) U (y) +V (x) V (y) , 
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The axial ratio of the ellipse is then given by the ratio of 
the Eigenvalues Ai,A2 determined by 



E - A F 
F G- A 







We find 



Ai , A2 



E + G± y / (E-G) 2 +4F 2 



(A8) 
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, - 1 ^ T ' . (AlO) 

Ai E + G+ y/(E- G) 2 +4F 2 

For the angle of any principal axis with the positive «-axis 
we find 

2F 

Although these equations get simpler if lensing properties 
are introduced (e. g. since = V( x j the last equation 
reads tan2y = 2u^/ (u^ — V( y ))) the complex notation 
is much more convenient. The equivalence of Eqs. (AlO, 
All) with Eqs. (12, 13) follows after a little computation if 
the definition of /j, Eq. (14) is fully expanded into normal 
derivatives 

("(*)-%)) +i ("(*) + "(y)) rA1 „, 
(«(*) + %)) + 1 (v (x) - u (y) ) 



References 

Ahlfors, L. V. 1966, Qusasiconformal Mappings, D. Van Nos- 

trand Company, Inc., Princeton, New Jersey. 
Breimer, T.G. 1994, A&A 283, 745 

Broadhurst, T.J., Taylor, A.N., Peacock, J. A. 1994, preprint 

no. 9406052 from astro-ph@babbage.sissa.it 
Carlberg, R.G., Yee, H.K.C., Ellington, E. 1994, preprint no. 

9402036 from astro-ph@babbage.sissa.it 
Carleson, L., Gamelin, T. W. 1993, Complex Dynamics, 

Springer- Verlag Heidelberg 
Fahlman, C.G., Kaiser, N., Squires, G., Woods, D. 1994, 

preprint no. 9402017 from astro-ph@babbage.sissa.it 
Fort, B., Mellier, Y. 1994, A&AR, 5, 239 

Imayoshi, Y., Taniguchi, M. 1992, An Introduction to Te- 

ichmiiller Spaces, Springer Verlag Tokyo 
Henrici, P. 1986 Applied and computational complex analysis, 

Vol III, §15.10, John Wiley and Sons, New York 
Kassiola, A., Kovner, I. 1993, ApJ, 417, 450 
Kaiser, N., Squires, G. 1993, ApJ, 404, 441 
Kneib, J.-P., Mathez, G., Fort, B., Mellier, Y., Soucail, G., 

Longaretti, P.-Y. 1994, A&A, 286, 701 
Miralda-Escude 1991, ApJ, 370,1 

Miralda-Escude, Babul, A. 1994, preprint no. 9405062 from 

astro-ph@babbage.sissa.it 
Schneider, P., Seitz, C. 1994, preprint no. 9407032 from astro- 

ph@babbage.sissa.it 
Schramm, T. 1994, A&A, 284,44 
Schramm, T., Kayser, R. 1987, A&A, 174, 361 
Schramm, T., Kayser, R. 1994, A&A in press 
Smail, I., Ellis, R.S., Fitchett, M.J. 1994 MNRAS, submitted, 

preprint, no. 9402048 from astro-ph@babbage.sissa.it 
Smail, I., Ellis, R.S., Fitchett, M.J. 1994 MNRAS, submitted, 

preprint, no. 9402049 from astro-ph@babbage.sissa.it 
Witt, H. J. 1990, A&A 236, 311 



This article was processed by the author using Springer- Verlag 
M? E X A&A style file L-AA version 3. 



-oc 


)oo< 


ooc 






00 O 

■0 c 

c 


1 07 OO < 

) 0< 


co c 

y \ ) V 

— ^ V 

Oo^ 


n c 

( 


7 O 

J OJ 

) 


■0 ( 

o <; 

c 

■oc 

-2 


) 0< 

) \Vo« 

} ^ O- 
)0O< 

-i 


30 . 

— O c 

^oo? c 
ooc 

-s 


' t 

? c 

DC 

1 


) 
1 
y 
"JO 

DO 

2 



X 



^ : -0 O O 



^ ^ ^ ^ 



0^ 

1 \o 



. 5- 



-0 . 5 



-1- 





^ ^ \ 



^ ^ ^ ^ 



^ ^ ^ ^ 



/ o o o ^ \ 
^ 
^ 





^ ^ ^ 




/ 



^000 



-1 . 5 



-1 



-0 . 5 



. 5 



1 . 5 



5- 



3 



2 



1- 



1 



